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ABSTRACT 



The Auxiliary Tracking System (ATS) is being developed to 
allow for an accurate underwater tracking capability that can 
be deployed anywhere in the world. A simulation was created 
to determine the impact upon the relative accuracy of the 
system due to a normally distributed noise. The simulation 
was used to look at the impact of range between targets and 
the depth spacing of the suspended transducers on the relative 
tracking accuracy. 
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I. 



INTRODUCTION 



Within the last several years, the Navy has seen a 
tightening of the purse strings from congress, and as a 
result, must find more cost effective ways to perform the 
various missions that are still required. One area that has 
long been a concern is the limited number of underwater 
tracking ranges. Alternatives to the development of another 
fixed range must be considered because of the financial 
impact. 

On the west coast, there are currently three areas where 
underwater tracking ranges are available. Washington, San 
Clemente, California, and Hawaii. With the limited number of 
sites and the other demands on range time for research and 
development, fleet ASW training opportunities are limited. 
Added to this is the fact that for many exercises, the units 
must transit to Hawaii simply because that is the only range 
available. The portable Auxiliary Tracking System (ATS) was 
conceived for these reasons. 

The ATS has been designed as a truly portable underwater 
tracking range. This thesis will look at the theoretical 
tracking accuracy of the system currently under engineering 
development and testing at the Naval Undersea Warfare 
Engineering Station (NUWES) . 
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II 



ATS CONCEPT 



A. GENERAL DESCRIPTION 

With the recent improvements made in electronics and the 
Global Positioning Satellite System (GPSS) , the dream of a 
portable tracking range now seems possible with current 
technology. The engineering development model consists of 
four underwater tracking buoys and a portable computer site. 
The buoys are deployed either from a ship or helicopter and 
are recoverable. For the four buoy system, deployment will be 
similar to that shown in Figure 1. The final system will 
consist of nine buoys to be deployed in various patterns based 
upon the needs of the system being tested. 




Figure 1. ATS Deployment Configuration 
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1. Buoy System 

The buoys used in the ATS will be free floating, self 
contained buoys. They will have a transducer suspended below 
at a predetermined depth (between 50 and 3 00 meters) used both 
for tracking the targets and locating the other units in the 
tracking range. The acoustic signal processor contained in 
each buoy will provide for the detection of the acoustic pings 
emitted by the other sources. The GPSS receiver in the buoy 
will provide the Latitude/Longitude and a very accurate time 
standard for the system. The on-board computer will control 
the GPSS receiver and acoustic signal processor, monitor the 
status of the buoy, perform a system up/down self test, time 
stamp and format the position data collected by the GPSS 
receiver and acoustic signal processor and send the data to 
the computer site via an RF data link. 

2. Computer Site 

The control station for the ATS range will be the 
computer site, which will be positioned on a ship located off 
the range, or on land if the range is set up in coastal 
waters. The computer site will be contained in a portable van 
allowing transport anywhere in the world either by ship or by 
cargo aircraft. The only support needed will be AC power, as 
the computer site will be self contained. The computer site 
will contain the RF data link receivers, RF voice 
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communications equipment, a GPSS receiver, the range tracking 
computers and peripherals, and other equipment as needed. 

B. SYSTEM OPERATION 

The concept embodied by the ATS is a "sing-around" self 
calibrating tracking range. This means that the buoys will 
locate each other. Each buoy has an acoustic pinger on the 
floating buoy and a transducer suspended below. In a two 
minute cycle, each buoy pinger will transmit an acoustic 
signal followed by each suspended transducer also transmitting 
an acoustic signal. The acoustic processor within each buoy 
will control the pinging cycle; a coded pulse will be used to 
identify the signal. The transit time of the signal will be 
known as each buoy will transmit a pulse at a specific time, 
and the reception time of the signal will be sent to the 
computer site via the RF link from the buoy that received the 
signal. From this, the transit time of the signal is known, 
and the range between the source and receiver is calculated by 
using the Sound Velocity Profile of the water mass. The 
signal from the buoy pinger to the transducer suspended below 
is used to calculate the depth of each transducer. The 
transit times from transducer to transducer are then used to 
calculate the slant ranges when the depth of each transducer 
has been determined. From this data, a range grid is set up 
as shown in Figure 2 . The location of each transducer and the 
range grid will be updated every two minutes by this method. 
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X 



Figure 2. ATS Range Grid 



with a relative grid established, the system now determines 
the positions of the targets on the range. The ATS has three 
distinct frequency bands to utilize, by incorporating a Shift 
Frequency, Shift Keying (SFSK) processor and transmitter. The 
high frequency will be used for the sing-around calibration 
already described, and the mid and low frequencies will be 
used to track the torpedo and the target. The SFSK processor 
and transmitter is similar to the system used on the Quinalt 
Shallow Water range at NUWES. 
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III. SIMULATION DESCRIPTION 



A. PROGRAM DESCRIPTION 

A simulation program was written to determine the effects 
of a normal error on the transit times received at the 
computer site to investigate the accuracy of the ATS. The 
effects of shadow zones were not considered in the simulation 
due to the relatively short distances between any source and 
receiver pair. The simulation program has three main 
sections; solving for transit times, solving for positions, 
and collecting statistics. 

1. Simulation Inputs 

The ATS simulation program requires three sets of 
inputs. The first set of inputs is the sound speed profile 
for the water mass. The program can accept up to 60 data 
points for the profile. The second set of inputs is the 
actual positions for the four buoys, the four suspended 
transducers and the two targets. The inputs are cartesian co- 
ordinate with buoy one being at the origin and buoy three 
establishing the X-Z plane. The third set of inputs is for 
the number of replications desired and the standard deviation 
of the noise for the systems. The standard deviation is a 
function of the distance between the source and receiver, and 
within the program is modeled by changing the travel times in 
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each of three groups; buoy to transducer, transducer to 
transducer, and target to transducer. 

2. Solving for Times 

Once the inputs have been entered, the program solves 
for the actual transit time between each source and receiver 
pair. This is done by taking the slant range between each 
source and receiver (Equation 3.1) and dividing by the 

Range=yJ ( Y^-Y^) ^ (3.1) 

effective sound speed between source and receiver. The 
effective sound speed is calculated by solving for the travel 
time for a vertically propagating wave in each layer. The 
travel times are summed between the source and receiver depths 
and the difference in depths is divided by the travel time to 
yield the effective sound speed between the source and 
receiver. The derivation of the equation is shown below and 
starts with the definition of the travel time (Equation 3.2). 

t=fdt=f—^ (3.2) 

J J c(z) 

So 

The ray path within a constant gradient is defined as the arc 
of a circle of radius equal to a (Equation 3.3 and 3.4) 

[Ref l:p. 402]. 

ds=odB (3.3) 
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(3.4) 



g COS0 

Substitution of Equation 3.4 into Equation 3.2 yields Equation 
3.5 and 3.6. The solution of the integral is shown in 
Equation 3.7 




d0 

COS0 



(3.5) 



t 



e 



1 

g 



j sec (0) dO 

0o 



(3.6) 



t=— [Inlsec (0) +tan (0) I] (3.7) 

9 

Separation of the terms yield Equation 3.8. Using snell's law 
[Ref 1: p. 401], as shown in Equation 3.9, and making the 
assumption that the sine terms cancel for a vertically 
propagating wave. Equation 3.8 can be reduced to Equation 
3.10. 



t= — ln| 
9 



l+sin0 

COS0 



cos0^ 

l+sin0^ 



(3.8) 



COS0q _ c(Zp) 
COS0 c(z) 



t= 




c{z^) 

{c(z^) 



(3.10) 
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the 



• By summing the travel times within each layer, 
propagation time is found. If within a layer, the gradient is 
equal to zero, the program takes the vertical distance within 
the layer and divides by the sound speed within the layer to 
determine the time t. The effective sound speed is now found 
(Equation 3.11). 

effspd=^^^ (3.11) 

Knowing the slant range and the effective sound speed between 
each source and receiver pair, the transit times for all pairs 
are determined (Equation 3.12). 



Tlme^^§rne 

effspd 



(3.12) 



3. Adding Noise 

When all of the transit times have been calculated and 
stored, the program begins the iteration cycle. The first 
step is the addition of the noise value to the transit times. 
This is done by taking a standard normal random variable and 
multiplying it by the transit time and the percentage value 
for the standard deviation. This results in a normal random 
variable with mean equal to zero and standard deviation equal 
to the given percentage of the transit time. This method was 
chosen because the transducer to transducer ranges and times 
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are typically 5 times greater than the pinger to suspended 
transducer times, and the error should be a function of range. 

4. Solving for Transducer Locations 

After the noise has been added to each of the transit 
times, the location of each transducer is calculated. The 
depth of each transducer is found first using the transit time 
from the buoy pinger to the transducer suspended below (i.e. 
buoy pinger one to transducer one) . From this transit time 
and using the length of the cable as a starting point, the 
depth is determined through an iterative process by 
calculating the effective sound velocity between the pinger 
and the estimated depth of the transducer. The solution is 
then found by multiplying the transit time by the effective 
sound velocity. The resulting depth is then compared with the 
estimated depth, and when the difference between the two 
values is less than one foot, the program continues, otherwise 
the resulting depth is used for the estimated depth and the 
calculation is repeated. Once the four transducer depths are 
established, the range coordinate system is established using 
transducer one as the origin of the X-Y grid and transducer 
three as the intersection with the horizontal axis (X-axis) . 
There are two transit times for each range calculation made 
(transit time from transducer one to three and from transducer 
three to one) , and the times are averaged. The slant range 
(SR) is then found by calculating the effective sound velocity 
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between the two depths and multiplying by the average transit 
time between the two transducers. The horizontal distance (H) 
is calculated using Equation 3.13. The position of transducer 

{SR)^-{Z^-Z2)^ ( 3 . 13 ) 

two is found by using the horizontal ranges from transducers 
one and three and using the GPSS position of buoy two to 
determine which of two unique solutions is correct (the 
intersection of two circles) . Transducer four is then located 
the same way using the positions and ranges from transducers 
two and three. 

5. Solving for Target Locations 

With the range coordinate system established and the 
location of all transducers found, we are positioned to 
calculate the location of the two targets. It is assumed that 
all four transducers receive the signal from each target. The 
location of the target is then calculated through a least 
squares method similar to one used by Larry Anderson at NUWES. 
The detailed procedures for the algorithm are presented below. 

Given the location of each transducer as (Xj, Yj, Z;) and 
the solution range between transducer i and the target (Rmi) 
from the transit time, the location of the Target (Xj, Yj, Zj) 
is calculated. is the measured slant range between target 
and transducer i, and R; is the calculated slant range between 
target and transducer i as given by Equation 3.14. The 
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Ri=^(X^-Xi)^+{Y^-Y^) 2+ (Zj,-Zj) 2 



(3.14) 



location of the target is then calculated with respect to 
minimizing the Sum of Squares of errors (SSE) , where the SSE 
is defined by Equation 3.15. The calculated slant ranges (Rj) 

SSE=Y, (3.15) 

i=l 



are linearized around the trial value Zj^) , shown in 

Equation 3.16, where | „ means evaluated at point "o”. The 



dR. 



dXn 



(Xr-X^ ) + 



dR, 



dY. 



(Y^-Yr ) + 



dR, 



dZ, 



(Z 3 .-Z™) (3.16) 



partial derivatives are listed in Equations 3.17 through 3.19. 





1 

II 


dXj. 


0 

1 

0 


dR, 


1 


dY^ 


0 

1 

0 


dR, 




dZj. 


0 

1 

0 



(3.17) 



(3.18) 



(3.19) 



The linearized model is substituted into Equation 3.15 and the 
SSE is minimized with respect to (X,., Yj, Zj) . The resulting 
normal equations are put in matrix form as shown in Equation 
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3.20. The individual matrices are shown in Equations 3.21 
through 3.23. Equation 3.20 is solved for the correction 
matrix (Equation 3.22) which when added to 



AX=B 



( 3 . 20 ) 



A = 



dR, 


2 




dR, 


dR, 


dX^ 


O 


E 


dY^ 


dXr 
o •' 



dR,. 



dR, 



dXj. 

dR, 



dR, 



dXn 



dYn 



dR, 



dR, 



dY. 



dZj 

dR, 



dR, 



dZn 



dR, 



dY. 



dR, 



dZn 



dZ, 



dXj. 



dR, 



dY. 



dR, 



dZ, 



( 3 . 21 ) 



X = 



X^-X^^ 

Z^-Z^ 



( 3 . 22 ) 



B = 






dR, 



dXj. 

dR, 



dY, 



( 3 . 23 ) 






yields the improved target position (X^., Y-j-, Z-r) . If the sum 
of the squares of the corrections is less than some value 
epsilon, Equation 3.24, the iteration is stopped. Otherwise, 



{x^-x^ ) 2+ ( y„-y„ ) 2+ (z.-z™ ) 2<e 



( 3 . 24 ) 



13 



the target position (X^, Y^, Zj) is used for the next iteration 
in place of (X^o/ Y^o, Zxo) and the procedure is repeated until 
50 iterations have been completed. 

6. Collecting Statistics 

When the positions of the targets have been solved 
for, the errors are calculated. The distance between the 
actual and solution location of each target (Equation 3.25), 
and the distance in the relative position of the targets with 
respect to each other (Equation 3.26) are calculated. 



Xoi is the X-value for the true position of the ith target, X,j 
is the X-value for the solved position of the ith target , TE; 
is the True error for the ith target, and RE is the relative 
error between the two targets. The values for TEj and RE are 
summed over the iteration loop and the programs final output 
is the mean and standard deviation of the true position error 
for each target and the relative position error between the 
two targets. 

B. SIMULATIONS CONDUCTED 

Three factors are considered as to the influence of the 
relative accuracy: Magnitude of the noise in the 




(3.25) 



RE=y/ ( • 



(3.26) 
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measurements, the proximity of target one to target two, and 
whether the transducers are at equal or staggered depths. 

1. Noise Values 

Noise values are entered into the simulation by the 
standard deviation being a percentage of the true value. 
Values examined range between 0.1 and 1.0 percent. 

2 . Target Locations 

The simulations were conducted for targets at close 
range (100 ft separation) and long range (3000 ft separation) 
to investigate the effects of target separation on the 
relative accuracy. 

3. Transducer Depths 

Some discussion has taken place as to the effects of 
staggering the transducer depths or having them all at an 
equal depth with the thought being would accuracy be gained by 
having the suspended transducers staggered at varying depths. 
Would this allow for better target resolution given the four 
slant ranges from the target to the transducers? Simulations 
were also conducted to draw a conclusion on the selection of 
transducers depths; equal depths or staggered depths with 
intervals of 100 feet. 
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IV. DATA ANALYSIS 



Each of the four following cases were run using the 
simulation program with the noise values ranging from 0 to 1 
percent. The performance goal for the ATS is 20 yards 
relative accuracy. The specifics of the four cases are as 
follows ; 

1. Case One: Transducers at equal depth (1000 feet) and 

targets at close range (100 feet) . 

2. Case Two: Transducers at equal depth and targets at 

long range (3000 feet) . 

3. Case Three: Transducers at staggered depths (100 foot 

intervals) and targets at close range. 

4. Case Four: Transducers at staggered depths and targets 

at long range. 

The results for the four cases are shown in Figure 3. The 
limit for the noise varies from 0.3 percent for cases two and 
four (targets at long range) to 0.4 percent for cases one and 
three (targets at close range) . 
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•- Case 1 Case 2 -*•- Case 3 Case 4 



Figure 3. Relative Error Plot 
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V. CONCLUSIONS AND RECOMMENDATIONS 

A. CONCLUSIONS 

The results from the simulation show that the system can 
not handle a lot of noise as currently designed. The 
staggering of the transducer depths did not appear to have an 
effect on the accuracy of the system. All four cases were 
relatively close for noise value where the relative accuracy 
exceeded 20 yards. 

B. RECOMMENDATIONS 

The accuracy assessment for the ATS is a very complex 
problem. This study has only focused on a few areas and a 
complete assessment of the system will require further study. 
From the initial results on the engineering development model, 
the outlook is very promising for the Navy to receive 
substantial benefits from its continued funding. Analysis of 
the data from the prototype buoys should reveal the magnitude 
of the noise, and the expected accuracy can then be 
determined. An additional area that should also be considered 
is to use telemetry data from the targets in the tracking 
problem. This could be utilized to provide the depth of the 
target and the velocity. By incorporating a Kalman Filter 
into the algorithm, it may be possible to reduce the impact 
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from a bad ping cycle. This would also necessitate looking at 
the accuracy of the depth sensors currently available. 

The Auxiliary Tracking System should continue to be 
developed. Potential benefits for the Navy warrant the 
investment now, as current ranges are taxed to keep up with 
the demands placed upon them by the many users. As new 
weapons systems are developed, the ATS will allow more 
realistic operational testing to be conducted by allowing the 
weapons systems to be evaluated in various ocean types. 
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APPENDIX A 



cc 

cc 

cc 

cc 

cc 



ATSSIM.FOR 



SIMULATION PROGRAM TO EVALUATE AUXILLARY TRACKING SYSTEM 
ABSOLUTE AND RELATIVE ACCURACY ASSESSMENTS BASED UPON A 
RANDOM NOISE WITH A GIVEN MEAN AND VARIANCE. 



REAL*8 SEED 

REAL*4 XBUOY, YBUOY, ZBUOY,XPHONE,YPHONE, ZPHONE,XTGT, YTGT, ZTGT, 

+ MEANl, VARl, MEAN2, VAR2, MEANR, VARR, SRANGE, RANGE, NOISE, 

+ TTIMEl, HTIMEl, BTIMEl, TTIME, HTIME, BTIME, A, B 
INTEGER*4 BUOYS, MAXSSV, NCOUNT 

DIMENSION XBUOY ( 4 ) , YBUOY ( 4 ) , ZBUOY ( 4 ) , XPHONE ( 4 ) , YPHONE ( 4 ) , 

+ ZPHONE(4), SRANGE(6) ,BTIME1(4,4) ,HTIME1(4,4) , 

+ PHONEX(4), PHONEY(4) , PHONEZ(4) , DEPTH ( 40) , SSPVEL( 40 ) , 

+ XTGT(2) ,YTGT(2) , ZTGT { 2 ) , TGTX ( 2 ) , TGTY ( 2 ) , TTIMEl ( 2 , 4 ) , 

+ TGTZ ( 2 ) , BTIME (4,4), HTIME (4,4), TTIME (2,4), 

+ TERROR(2) ,TESQ(2) ,TESUM(2) ,TRANGE(2,4) 

MAXSSV = 40 
BUOYS = 4 
ZCABLE = 1000.0 

CALL INPUT (NCOUNT, BTN, STN, TTN, SEED) 

CALL PROFILE (MAXSSV, NUMSSV, DEPTH, SSPVEL) 

CALL POSITS (XBUOY, YBUOY, ZBUOY, XPHONE, YPHONE, ZPHONE, 

+ XTGT, YTGT, ZTGT) 

CC CALCULATE TRANSIT TIMES FOR BUOYS 
DO 20 I = 1, 4 
DO 10 J = 1, 4 

RANGE = SQRT( (XBUOY(I)-XPHONE(J) )**2+(YBUOY(I)-YPHONE(J) ) **2 + 
+ (ZBUOY(I) - ZPHONE(J) )**2) 

CALL HMVELOC(ZBUOY(I) ,ZPHONE(J) , NUMSSV, DEPTH, SSPVEL, HMVEL ) 
BTIME ( I, J) = RANGE/HMVEL 
10 CONTINUE 
20 CONTINUE 

DO 50 I = 1, 4 
DO 40 J = 1, 4 

IF (I .EQ. J) GOTO 30 
RANGE = SQRT( (XPHONE(I)-XPHONE(J) )**2+ 

+ (YPHONE(I)-YPHONE( J) ) **2+ ( ZPHONE ( I ) -ZPHONE ( J ) )**2) 

CALL HMVELOC( ZPHONE (I) , ZPHONE(J), NUMSSV, DEPTH, SSPVEL, HMVEL) 
HTIME(I,J) = RANGE/HMVEL 
30 CONTINUE 

40 CONTINUE 
50 CONTINUE 

DO 70 I = 1, 2 
DO 60 J = 1, 4 

RANGE = SQRT( (XTGT(I)-XPHONE(J) )**2+(YTGT(I)-YPHONE(J) )**2+ 

+ (ZTGT(I) - ZPHONE (J) ) **2) 

CALL HMVELOC(ZTGT(I) , ZPHONE (J) , NUMSSV, DEPTH, SSPVEL, HMVEL ) 

TTIME ( I, J) = RANGE/HMVEL 
60 CONTINUE 
70 CONTINUE 
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DO 900 KOUNT=l, NCOUNT 
79 CONTINUE 
CC ADD NOISE 

DO 71 I = 1, 2 
J = I + 2 

CALL URAND{A,B,SEED) 

BTIME1(I,I) = BTIME(I,I) + ( A*BTN*BTIME ( I , I ) ) 
BTIME1(J,J) = BTIME(J,J) + (B*BTN*BTIME( J, J) ) 

71 CONTINUE 

DO 73 I = 1, 4 
DO 72 J = 1, 2 
K = J + 2 

CALL URAND (A, B, SEED) 

HTIME1(I,J) = HTIME(I,J) + { A*STN*HTIME { I , J) ) 
HTIME1{I,K) = HTIME{I,K) + { B*STN*HTIME ( I , K) ) 

72 CONTINUE 

73 CONTINUE 

DO 75 J = 1, 4 
CALL URAND (A, B, SEED) 

TTIME1(1,J) = TTIME(1,J) + { A*TTN*TTIME ( 1 , J ) ) 
TTIME1(2,J) = TTIME(2,J) + ( B*BTN*TTIME { 2 , J) ) 



75 CONTINUE 

CC CALC BUOY POSITS 
DO 80 I = 1, 4 
PHONEX(I) = XBUOY(I) 

PHONEY (I) = YBUOY(I) 

80 CONTINUE 

CALL HYDEPTH (ZCABLE, ZBUOY, NUMSSV, DEPTH, 

+ SSPVEL, HMVEL, BTIMEl, PHONEZ) 

CALL RPHONE (SRANGE, HTIMEl, PHONEZ, NUMSSV, DEPTH, SSPVEL, HMVEL) 
CALL TWOD (SRANGE, PHONEX, PHONEY, PHONEZ, YBUOY) 

CC CALC TARGET POSITS 

DO 100 J = 1, 2 
TGTZ(J) = 400 
TGTX(J) = 2000 
TGTY(J) = 2000 

CALL TGTRANGE (TTIME1,TGTZ, ZPHONE, NUMSSV, DEPTH, SSPVEL, HMVEL, 

+ TRANGE,J) 

CALL THREED (J, TGTX, TGTY, TGTZ, PHONEX, PHONEY, PHONEZ ,TRANGE) 

100 CONTINUE 
CC CALC TRUE ERRORS 

ERFIX = 0.0 
DO 150 J = 1,2 

TERROR(J) = SQRT( (XTGT(J)-TGTX(J) ) **2+(YTGT(J)-TGTY(J) ) **2+ 

+ (ZTGT(J)-TGTZ(J) )**2) 

IF (TERROR(J) .GT. 1000.0) ERFIX = 1.0 
150 CONTINUE 

IF (ERFIX .EQ. 1) THEN 
GOTO 79 
ELSE 

DO 200 J =1,2 

TESQ(J) = TESQ(J) + TERROR(J)**2 
TESUM(J) = TESUM(J) + TERROR(J) 

200 CONTINUE 

CC CALC RELATIVE ERRORS 

RERROR = SQRT( (XTGT(1)-TGTX(1)-XTGT(2)+TGTX(2) )**2+ 

+ (YTGT(1)-TGTY(1)-YTGT(2)+TGTY(2) )**2+ 

+ (ZTGT(1)-TGTZ(1)-ZTGT(2)+TGTZ(2) ) **2) 

RESQ = RESQ + RERROR* *2 
RESUM = RESUM + RERROR 
ENDIF 
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900 CONTINUE 

CC CALC STATS ON TRUE ERRORS 
CC TGTl RESULTS 

MEANl = TESUM(l)/NCOUNT 

VARl = SQRT(ABS(NCOUNT*TESQ(l)-TESUM(l)**2) )/(NCOUNT*(NCOUNT-l) ) 
CC TGT2 RESULTS 

MEAN2 = TESUM(2)/NCOUNT 

VAR2 = SQRT(ABS(NCOUNT*TESQ(2)-TESUM{2)**2) )/{NCOUNT*(NCOUNT-l) ) 
CC CALC STATS ON RELATIVE ERRORS 
CC RELATIVE RESULTS 

MEANR = RESUM/NCOUNT 

VARR = SQRT(ABS(NCOUNT*RESQ-RESUM**2))/(NCOUNT*(NCOUNT-l) ) 

CALL OUTPUT (MEANl, VARl, MEAN 2 , VAR2, MEANR, VARR, 

+ BTN, STN,TTN,NCOUNT) 

END 

* INPUT SUBROUTINE TO ENTER IN SIMULATION VALUES * 

SUBROUTINE INPUT (NCOUNT, BTN, STN, TTN, SEED) 

REAL* 8 SEED 

OPEN (UNIT = 1, FILE = 'INPUT.DAT') 

READ (1,*) NCOUNT, BTN, STN, TTN, SEED 
CLOSE (UNIT = 1, STATUS = 'KEEP') 

RETURN 

END 

* PROFILE SUBROUTINE TO ENTER SOUND VELOCITY PROFILE * 

**************************************************************** 

SUBROUTINE PROFILE (MAXSSV, NUMSSV, DEPTH, SSPVEL) 

REAL* 4 DEPTH, SSPVEL 
INTEGER*4 MAXSSV, NUMSSV 
DIMENSION DEPTH(40), SSPVEL(40) 

OPEN (UNIT = 1, FILE = 'SSPROF.DAT') 

READ (1,*) NUMSSV 
IF (NUMSSV .GT. MAXSSV) THEN 

WRITE(*,*) ' MAXIMUM NUMBER OF SSP DATA POINTS EXCEEDED' 
ENDIF 

DO 10 I = 1, NUMSSV 

READ (1,*) DEPTH(I), SSPVEL(I) 

10 CONTINUE 

CLOSE (UNIT = 1, STATUS='KEEP' ) 

RETURN 

END 



**************************************************************** 

* POSITS SUBROUTINE TO ENTER BUOY, TRANSDUCER AND TARGET * 

**************************************************************** 



SUBROUTINE POSITS ( XBUOY , YBUOY , ZBUOY , XPHONE , YPHONE , ZPHONE , XTGT, 
+ YTGT,ZTGT) 

REAL* 4 XBUOY, YBUOY, ZBUOY, XPHONE, YPHONE, ZPHONE 
DIMENSION XBUOY ( 4 ) , YBUOY ( 4 ) , ZBUOY ( 4 ) , XPHONE ( 4 ) , YPHONE ( 4 ) , 

+ ZPHONE (4) ,XTGT(2) ,YTGT(2) ,ZTGT(2) 

OPEN (UNIT = 2, FILE = 'BUOYLOC.DAT') 

DO 10 I = 1, 4 

READ ( 2 , * ) XBUOY ( I ) , YBUOY ( I ) , ZBUOY ( I ) , XPHONE ( I ) , YPHONE ( I ) , 

+ ZPHONE(I) 



10 CONTINUE 

OPEN (UNIT = 1, FILE = 'TARGET.DAT') 
DO 20 J = 1, 2 
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ooo oooooo o o oo 



READ(1,*) XTGT(J), YTGT(J), ZTGT(J) 

20 CONTINUE 
RETURN 
END 

**************************************************************** 

* HMVELOC SUBROUTINE TO CALCULATE MEAN SOUND VELOCITY * 

SUBROUTINE HMVELOC (ZPING, ZRCVR, NUMSSV, DEPTH, SSPVEL, HMVEL) 
REAL*4 ZPING, ZRCVR, DEPTH, SSPVEL, HMVEL, GRADNT, HMVl, HMV2,T 
INTEGER*4 I, J, K, INIT, FINAL, NUMSSV 
DIMENSION DEPTH ( NUMSSV) , SSPVEL (NUMSSV) 

ORDER THE DEPTH OF THE RECEIVER AND HYDROPHONE SO THAT THE 
RECEIVER IS THE DEEPER OF THE TWO DEPTHS. 

T = 0.0 
INTE = 0 

IF (ZPING .GT. ZRCVR) THEN 
TEMP = ZPING 
ZPING = ZRCVR 
ZRCVR = TEMP 
INTE = 1 
ENDIF 

INTERPOLATE THE SOUND SPEED AT THE PINGER DEPTH 
DO 1 I = 2, NUMSSV 

IF (ZPING .GT. DEPTH(I)) GOTO 1 

GRADNT = (SSPVEL(I) - SSPVEL ( I-l ))/ (DEPTH ( I ) - DEPTH (I-l)) 
HMVl = SSPVEL (I-l) + GRADNT* (ZPING-DEPTH( I-l ) ) 

INIT= I 
GO TO 2 

1 CONTINUE 

INTERPOLATE THE SOUND SPEED AT THE HYDROPHONE DEPTH • 

2 DO 3 J = 2, NUMSSV 
IF (ZRCVR .GT. DEPTH (J)) GOTO 3 

GRADNT = (SSPVEL(J) - SSPVEL (J-1 ))/ (DEPTH(J) - DEPTH (J-1)) 
HMV2 = SSPVEL (J-1) + GRADNT* (ZRCVR-DEPTH( J-1) ) 

FINAL = J-1 
GO TO 4 

3 CONTINUE 

CALCULATE THE HARMONIC SOUND SPEED BETWEEN THE PINGER AND THE 
HYDROPHONE 

4 IF (FINAL .LT. INIT) THEN 
SOURCE AND RECEIVER IN SAME GRADIENT 

TAKE AVERAGE OF TWO SOUND SPEEDS OF SOURCE AND RECEIVER 

GRADNT = (HMV2 - HMVl) / (ZRCVR - ZPING) 

IF (GRADNT .NE. 0.0) T = LOG (HMV2 /HMVl ) /GRADNT 

IF (GRADNT .EQ. 0.0) T = (ZRCVR - ZPING)/HMV2 

HMVEL = T/(ZRCVR-ZPING) 

ELSEIF (FINAL .EQ. INIT) THEN 

SOURCE AND RECEIVER SPERATED BY ONE DEPTH GRADIENT 

GRADNT = (SSPVEL(INIT) - HMVl )/ (DEPTH ( INIT) - ZPING) 

IF (GRADNT .NE. 0.0) T = LOG ( SSPVEL ( INIT) /HMVl ) /GRADNT 

IF (GRADNT .EQ. 0.0) T = (DEPTH( INIT) -ZPING) /SSPVEL( INIT) 

GRADNT = (HMV2-SSPVEL(FINAL) )/ (ZRCVR - DEPTH ( FINAL ) ) 

IF (GRADNT .NE. 0.0) T = T + LOG(HMV2/SSPVEL(FINAL) ) /GRADNT 

IF (GRADNT .EQ. 0.0) T = T + ( ZRCVR-DEPTH ( FINAL) ) /HMV2 

HMVEL = ABS( (ZRCVR-ZPING)/T) 

ELSE 
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. c 

C SOURCE AND RECEIVER SEPERATED BV MORE THAN ONE DEPTH 

C 

GRADNT = (SSPVEL(INIT) - HMVl )/( DEPTH ( INIT) - ZPING) 

IF (GRADNT .NE. 0.0) T = LOG ( SSPVEL ( INIT) /HMVl ) /GRADNT 
IF (GRADNT .EQ. 0.0) T = ( DEPTH ( INIT ) -ZPING ) /SSPVEL ( INIT ) 

DO 5 K= INIT + 1, FINAL 

GRADNT = (SSPVEL(K)-SSPVEL(K-l) )/(DEPTH(K)-DEPTH(K-l) ) 

IF (GRADNT .NE. 0.0) T = T+LOG ( SSPVEL ( K) /SSPVEL ( K-1 )) /GRADNT 
IF (GRADNT .EQ. 0.0) T = T+ ( DEPTH ( K) -DEPTH (K-1 )) /SSPVEL (K ) 

5 CONTINUE 

GRADNT = (HMV2-SSPVEL ( FINAL) )/(ZRCVR - DEPTH ( FINAL ) ) 

IF (GRADNT .NE. 0.0) T = T + LOG (HMV2 /SSPVEL ( FINAL) ) /GRADNT 
IF (GRADNT .EQ. 0.0) T = T + ( ZRCVR-DEPTH ( FINAL )) /HMV2 
HMVEL = ABS( (ZRCVR-ZPING)/T) 

ENDIF 

IF (INTE .EQ. 1) THEN 
TEMP = ZPING 
ZPING = ZRCVR 
ZRCVR = TEMP 
ENDIF 
RETURN 
END 



* HYDEPTH SUBROUTINE TO CALCULATE TRANSDUCER DEPTHS * 

SUBROUTINE HYDEPTH (ZCABLE, ZBUOY, NUMSSV, DEPTH, 

+ SSPVEL, HMVEL, BTIME, PHONEZ) 

CC SUBROUTINE TO CALCULATE THE DEPTH OF THE HYDROPHONE BELOW THE 
CC BUOY • 

REAL*4 ZCABLE, ZBUOY, DEPTH, SSPVEL, HMVEL, BTIME , ZOLD , PHONEZ 
INTEGER*4 NUMSSV 

DIMENSION DEPTH(40), SSPVEL ( 40 ), PHONEZ ( 4 ) , BTIME ( 4 , 4 ) , ZBUOY ( 4 ) 
DO 5 I = 1, 4 
ZB = ZBUOY(I) 

ZOLD = ZCABLE 
DO 3 KOUNT = 1, 10 

CALL HMVELOC(ZB, ZOLD, NUMSSV, DEPTH, SSPVEL, HMVEL) 

ZNEW = HMVEL * BTIME (I, I) + ZB 

IF (ABS(ZNEW-ZOLD) .LT. 1.0) GO TO 4 
ZOLD = ZNEW 

3 CONTINUE 

4 PHONEZ (I) = ZNEW 

5 CONTINUE 
RETURN 
END 



****★***★★*******************★*★*★★★*★★★**★*★**★*★★*★*★★★*★*★*★* 

* RPHONE SUBROUTINE TO CALCULATE RANGES BETWEEN TRANSDUCERS* 

★**★*★★★**★★★*★*★**★★★**★★★*******★★***★**★**★***★*★★★★*★*★****★ 

SUBROUTINE RPHONE (RANGE, HTIME, PHONEZ, NUMSSV, DEPTH, SSPVEL, 
+ HMVEL) 



CC 

CC 


RANGE(l) = PHONE 


1 


TO 


PHONE 


2 


CC 


RANGE (2) = PHONE 


3 


TO 


PHONE 


2 


CC 


RANGE (3) = PHONE 


1 


TO 


PHONE 


4 


CC 


RANGE (4) = PHONE 


3 


TO 


PHONE 


4 


CC 


RANGE(5) = PHONE 


2 


TO 


PHONE 


4 
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,HMVEL,T 
DEPTH (40) 



CC RANGE (6) = PHONE 1 TO PHONE 3 

REAL*4 RANGE, HTIME, PHONEZ, DEPTH, SSPVEL,! 

INTEGER NUMSSV 

DIMENSION RANGE (6) , HTIME(4,4), PHONEZ(4), 

+ SSPVEL(40) 

DO 1 I = 1, 4 
1 CONTINUE 

T = (HTIME (2,1) - HTIME(1,2)/3.0)*1.5 

CALL HMVELOC(PHONEZ(l) , PHONEZ(2), NUMSSV, DEPTH, SSPVEL, 
+ HMVEL) 

RANGE (1) = T * HMVEL 
T = 0.5*HTIME(2,3) + 0. 5*HTIME ( 3 , 2 ) 

CALL HMVELOC(PHONEZ(2) , PHONEZ(3), NUMSSV, 

+ HMVEL) 

RANGE (2) = T * HMVEL 
T = 0.5*HTIME(1,4) + 0 . 5*HTIME ( 4 , 1 ) 

CALL HMVE LOC ( PHONE Z ( 1 ) , PHONE Z(4), NUMSSV, 

+ HMVEL) 

RANGE (3) = T * HMVEL 

T = (HTIME(3,4) - HTIME ( 4 , 3 ) /3 . 0 ) *1 . 5 
CALL HMVELOC (PHONEZ (3), PHONEZ (4), NUMSSV, 

+ HMVEL) 

RANGE(4) = T * HMVEL 

T = 0.75*HTIME(2,4) + 0 . 25*HTIME ( 4, 2 ) 

CALL HMVELOC(PHONEZ(2), PHONEZ(4), NUMSSV, 

+ HMVEL) 

RANGE(5) = T * HMVEL 

T = 0.25*HTIME(1,3) + 0 . 75*HTIME ( 3 , 1 ) 

CALL HMVELOC(PHONEZ(l) , PHONEZ (3), NUMSSV, 

+ HMVEL) 

RANGE(6) = T * HMVEL 
RETURN 
END 



DEPTH, SSPVEL, 



DEPTH, SSPVEL, 



DEPTH, SSPVEL, 



DEPTH, SSPVEL, 



DEPTH, SSPVEL, 



* * * * * Him * * * * * UK * * * ******** K * ****** * * * UK * * * * * Him * * * *11 * * * * * UK * UK'* * 



TWOD SUBROUTINE TO CALCULATE TRANSDUCER POSITIONS 



**************************************************************** 



SUBROUTINE TWOD (RANGE, XP, YP, ZP, BUOYY) 

REAL*4 RANGE, BUOYY 

REAL*4 XP, YP, ZP, DX, DY, ANS, XSQRl, XSQR3, 

+ YSQRl, YSQR3, RSQRl, RSQR3, RHS, D I SCRIM, THETA, PHI 
DIMENSION BUOYY(4), XP(4), YP(4), ZP(4), RANGE(6), 

+ ANS(2,2) 

XP(3)=SQRT(RANGE(6)**2 - ( ZP( 1 ) -ZP ( 3 ) ) **2 ) 

YP(3)=0.0 

DX=XP(1)-XP(3) 

DY=YP(1)-YP(3) 

RSQR1=RANGE ( 1 ) **2- ( ZP ( 1 ) -ZP ( 2 ) ) * *2 

RSQR3=RANGE ( 2 ) * *2- ( ZP ( 2 ) -ZP ( 3 ) ) * *2 

XSQR1=XP(1)**2 

XSQR3=XP(3)**2 

YSQR1=YP(1)**2 

YSQR3=YP(3)**2 

RHS = (XSQRl + YSQRl - XSQR3 - YSQR3 +RSQR3 - RSQRl) *0.5 
IF (DY .EQ. 0.0) THEN 
THETA = RHS/DX 
A = 1 

B = -2 * YP(1) 

C = (THETA)**2 + XSQRl + YSQRl - (2 * THETA * XP ( 1 ) ) - RSQRl 
DISCRIM = B**2 - (4 * A * C) 
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(-B + SQRT(ABS(DISCRIM) ) ) / (2*A) 
(-B/A - ANS(1,2) ) 

THETA 

THETA 



ANS (1,2) 

ANS(2,2) 

ANS (1,1) 

ANS (2,1) 

ELSE 

THETA = RHS/DY 
PHI = DX/DY 
A = 1.0 + PHI**2 

B = 2* (YP(1)*PHI-THETA*PHI-XP(1) ) 

C = XSQRl + (THETA)**2 + YSQRl - RSQRl - 2.0 * THETA * YP(1) 



B**2 - (4 * A * C) 

(-B +SQRT(DISCRIM) )/(2*A) 
(-B/A - ANS(1,1) ) 

- PHI*ANS(1,1) 

- PHI*ANS(2,1) 



THETA 

THETA 



- BUOYY(2)) .LT. 1000) THEN 



ANS (2,1) 
ANS (2,2) 



DISCRIM 
ANS(1, 1) 

ANS (2,1) 

ANS (1,2) 

ANS(2,2) 

ENDIF 

IF (ABS(ANS(1,2) 

XP(2) = ANS (1,1) 

YP(2) = ANS (1,2) 

ELSE 
XP(2) 

YP(2) 

ENDIF 
DX=XP(2)-XP(3) 

DY=YP(2)-YP(3) 

RSQR1=RANGE ( 5 ) ** 2 - ( ZP ( 2 ) -ZP ( 4 ) ) **2 
RSQR3=RANGE ( 4 ) **2- ( ZP ( 3 ) -ZP ( 4 ) ) **2 
XSQR1=XP(2)**2 
XSQR3=XP(3)**2 
YSQR1=YP(2)**2 
YSQR3=YP(3)**2 
RHS = (XSQRl + YSQRl - 
IF (DY .EQ. 0.0) THEN 
THETA = RHS/DX 
A = 1 

B = -2 * YP(2) 

C = (THETA) **2 + 

DISCRIM = B**2 - 
(-B + 

(-B/A 
THETA 
THETA 



XSQR3 - YSQR3 +RSQR3 - RSQRl) *0.5 



XSQRl + YSQRl - (2 * THETA * XP(2)) - RSQRl 
(4 * A * C) 

SQRT(DISCRIM) ) / (2*A) 

- ANS(1,2)) 



ANS (1,2) 

ANS (2, 2) 

ANS (1,1) 

ANS (2,1) 

ELSE 

THETA = RHS/DY 
PHI = DX/DY 
A = 1.0 + PHI**2 

B = 2* (YP(2) *PHI-THETA*PHI-XP(2) ) 

C = XSQRl + (THETA)* *2 + YSQRl - RSQRl - 2.0 * THETA * YP(2) 



THETA 

THETA 



DISCRIM 
ANS (1,1) 

ANS (2,1) 

ANS (1,2) 

ANS (2,2) 

ENDIF 

IF (ABS(ANS(1,2) - 
XP(4) = ANS (1,1) 
YP(4) = ANS (1,2) 

ELSE 



B**2 - (4 * A * C) 

(-B +SQRT(ABS(DISCRIM) ) )/ (2*A) 
(-B/A - ANS (1,1) ) 

- PHI*ANS(1,1) 

- PHI*ANS(2,1) 



BUOYY(4)) .LT. 1000) THEN 
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XP(4) = ANS(2,1) 

YP(4) = ANS(2,2) 

ENDIF 

RETURN 

END 

★★************************************************************** 

* TGTRANGE SUBROUTINE TO CALCULATE RANGE TO TARGETS * 

**************************************************************** 

SUBROUTINE TGTRANGE (TTIME,TGTZ, ZPHONE, NUMSSV, DEPTH , SSPVEL , 

+ HMVEL, TRANGE, J) 

REAL* 4 TGTZ , Z PHONE , DEPTH , SSPVEL , HMVEL , TRANGE , TTIME 
INTEGER* 4 J, NUMSSV 

DIMENSION TGTZ (2) , ZPHONE ( 4 ) ,DEPTH ( 40) , SSPVEL ( 40) , TRANGE ( 2 , 4 ) , 

+ TTIME (2, 4) 

DO 10 I = 1, 4 

CALL HMVELOC(TGTZ( J) , ZPHONE ( I ) , NUMSSV, DEPTH, SSPVEL, HMVEL) 
TRANGE (J, I) = TTIME (J, I) * HMVEL 
10 CONTINUE 
RETURN 
END 

**************************************************************** 

* THREED SUBROUTINE TO CALCULATE TARGET POSITIONS * 

* ★★★★*★***★★★*★★★*★★★*★*★★★★★★★*★★★★***★***★*★***★*★★★★*★*★★★★★★ 

SUBROUTINE THREED ( J, TGTX, TGTY, TGTZ , PHONEX, PHONEY , PHONE Z , RM) 
REAL* 4 TGTX , TGTY , TGTZ , PHONEX , PHONEY , PHONEZ , PRPX , PRPY , PRPZ , 

+ XTO, YTO, ZTO,Al, A2,A3,A4,A5,A6,A7,A8,A9,RO,DR,RM 
INTEGER*4 J, IN 

DIMENSION TGTX(2) , TGTY (2 ) , TGTZ (2 ) , PHONEX ( 4 ) , PHONEY ( 4 ) , 

+ PHONEZ(4) ,RO(4) ,RM(2,4) 

IN = 0 
5 CONTINUE 
XTO=TGTX(J) 

YTO=TGTY( J) 

ZTO=TGTZ( J) 

A1=0.0 

A2=0.0 

A3=0.0 

A4=0.0 

A5=0.0 

A6=0.0 

A7=0.0 

A8=0.0 

A9=0.0 



DO 10 


1 = 


1, 4 




RO(I) 


— 


SQRT( (XTO-PHONEX(I) ) **2+ (YTO-PHONEY ( I ) ) **2+ 
(ZTO-PHONEZ(I) )**2) 


DR = 


RM(J,I) - RO(I) 




PRPX 


= 


(XTO - PHONEX (I)) 


/ RO(i) 


PRPY 


= 


(YTO - PHONEY (I)) 


/ RO(i) 


PRPZ 


= 


(ZTO - PHONEZ(I)) 


/ RO(i) 


A1 = 


A1 


+ PRPX* PRPX 




A2 = 


A2 


+ PRPX*PRPY 




A3 = 


A3 


+ PRPX*PRPZ 




A4 = 


A4 


+ PRPY* PRPY 




A5 = 


A5 


+ PRPY* PRPZ 




A6 = 


A6 


+ PRPZ*PRPZ 




A7 = 


A7 


+ DR* PRPX 




AS = 


AS 


+ DR*PRPY 
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A9 = A9 + DR*PRPZ 
10 CONTINUE 

DETAO = Al* (A4*A6-A5*A5)-A2* (A2*A6-A5*A3)+A3* (A2*A5-A4*A3) 
DETAl = A7*(A4*A6-A5*A5)-A8*(A2*A6-A5*A3)+A9*(A2’‘A5-A4*A3) 
DETA2 = Al*(A8*A6-A9*A5)-A2*(A7*A6-A9*A3)+A3*(A7*A5-A8*A3) 
DETA3 = Al* (A4*A9-A5*A8)-A2* (A2*A9-A5*A7)+A3* (A2*A8-A4*A7) 
TGTX(J) = XTO + DETAl/DETAO 

TGTY(J) = YTO + DETA2 /DETAO 

TGTZ(J) = ZTO + DETA3 /DETAO 

IF (( (DETA1**2+DETA2**2+DETA3**2) /DETAO) .LT. 1.0 ) THEN 
GO TO 900 
END IF 

IN = IN + 1 
IF (IN .GT. 25) THEN 
GO TO 900 
ENDIF 
GOTO 5 
900 RETURN 
END 



* OUTPUT SUBROUTINE TO OUTPUT DATA RESULTS * 

SUBROUTINE OUTPUT (MEANl , VARl , MEAN2 , VAR2 , MEANR, VARR, 

+ BTN,STN,TTN,NCOUNT) 

REAL *4 MEANl , VARl , MEAN2 , VAR2 , MEANR, VARR, BTN, STN, TTN 
INTEGER*4 NCOUNT 

OPEN (UNIT = 8, FILE = 'OUTPUT3.DAT') 



WRITE(8, *) 
WRITE (8, 1000) 
WRITE (8, 1001) 
WRITE (8, 1002) 
WRITE (8,*) ' 

WRITE (8, 1003) 

1000 FORMAT 

1001 FORMAT 

1002 FORMAT 

1003 FORMAT 
RETURN 
END 



MEAN 



STAND DEV 



MEANl, VARl 
MEAN 2, VAR2 
MEANR, VARR 
NUMBER OF ITERATIONS, 
NCOUNT, BTN, STN, TTN 
TRUE ERROR FOR TARGET 1 = 
TRUE ERROR FOR TARGET 2 = 
RELATIVE ERROR 



( 

( 

( 

( 3X,I6,15X, F8.4, F8.4, F8.4) 



BUOY 

' ,F8.3, 
' ,F8.3, 
' ,F8.3, 



TRANS 

F8.4) 

F8.4) 

F8.4) 



TARGET 



•k-k-k-kicic-kiticicic-kit-kic-kit-kit-kicTkicicicitic-k-k-kiciticicicicit-kiciciciticiticicic-kiciciiiciciciciticiciciciticicic 

* URAND SUBROUTINE TO GENERATE RANDOM NUMBERS * 

**************************************************************** 

SUBROUTINE URAND (A, B, SEED ) 

REAL*8 SEED, PI, Rl, R2 
PI=3. 14159265358979 
R1=(SEED+PI) **504.D-2 
R1=R1-DINT(R1) 

SEED=R1 

R2=(SEED+PI) **504.D-2 
R2=R2-DINT(R2) 

SEED=R2 

R=SQRT(-2*LOG(Rl) ) 

T=2*PI*R2 

A=R*SIN(T) 

B=R*COS(T) 

RETURN 

END 
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APPENDIX B 



The following is the output data from the simulations 
conducted for each of the four cases discussed in Chapter IV. 



Case number 1: 



BUOYLOC . DAT 



0.0, 


0.0, 


5.0, 


0.0, 


0.0, 


1000.0 


3000.0, 


5196.0, 


5.0, 


3000.0, 


5196.0, 


1000.0 


6000.0, 


0.0, 


5.0, 


6000.0, 


0.0, 


1000.0 


9000.0, 


5196.0, 


5.0, 


9000.0, 


5196.0, 


1000.0 



TARGET . DAT 

2036.0, 2036.0, 400.0, 

1965.0, 1965.0, 400.0, 

OUTPUT . DAT 

TRUE ERROR FOR TARGET 1 
TRUE ERROR FOR TARGET 2 
RELATIVE ERROR 

NUMBER OF ITERATIONS, 
1000 



MEAN STAND DEV 

= .002 .0000 

= .001 .0000 

= .001 .0000 

BUOY TRANS TARGET 

.0000 .0000 .0000 



TRUE ERROR FOR TARGET 1 
TRUE ERROR FOR TARGET 2 
RELATIVE ERROR 

NUMBER OF ITERATIONS, 
1000 



MEAN 

19.829 

19.695 

14.982 



STAND DEV 
. 0122 
. 0122 
.0084 



BUOY 

.0010 



TRANS 
. 0010 



TARGET 
, 0010 



TRUE ERROR FOR TARGET 1 
TRUE ERROR FOR TARGET 2 
RELATIVE ERROR 

NUMBER OF ITERATIONS, 
1000 



MEAN 

39.791 

39.564 

30.101 



STAND DEV 
.0247 
.0248 
.0169 



BUOY 

.0020 



TRANS 

.0020 



TARGET 

.0020 



TRUE ERROR FOR TARGET 1 
TRUE ERROR FOR TARGET 2 
RELATIVE ERROR 

NUMBER OF ITERATIONS, 
1000 



MEAN 



STAND DEV 



= 60.054 .0380 

= 59.780 .0383 

= 45.545 .0260 

BUOY TRANS TARGET 
.0030 .0030 .0030 
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TRUE ERROR FOR TARGET 1 
TRUE ERROR FOR TARGET 2 
RELATIVE ERROR 

NUMBER OF ITERATIONS, 
1000 



MEAN STAND DEV 
= 80.887 .0530 

= 80.616 .0537 

= 61.654 .0362 

BUOY TRANS TARGET 
.0040 .0040 .0040 



TRUE ERROR FOR TARGET 1 
TRUE ERROR FOR TARGET 2 
RELATIVE ERROR 

NUMBER OF ITERATIONS, 
1000 



MEAN 

= 102.065 
= 102.172 
= 78.586 



STAND DEV 
.0688 
.0718 
.0486 



BUOY 

.0050 



TRANS 
. 0050 



TARGET 

.0050 



TRUE ERROR FOR TARGET 1 
TRUE ERROR FOR TARGET 2 
RELATIVE ERROR 

NUMBER OF ITERATIONS, 
1000 



MEAN 

= 123.875 
= 123.837 
= 98.655 



STAND DEV 
.0910 
.0939 
.0726 



BUOY 

.0060 



TRANS 

.0060 



TARGET 

.0060 



TRUE ERROR FOR TARGET 1 
TRUE ERROR FOR TARGET 2 
RELATIVE ERROR 

NUMBER OF ITERATIONS, 
1000 



MEAN 

= 145.096 

= 146.658 

= 118.392 



STAND DEV 
.1044 
. 1133 
.0931 



BUOY 
. 0070 



TRANS 
. 0070 



TARGET 
. 0070 



TRUE ERROR FOR TARGET 1 
TRUE ERROR FOR TARGET 2 
RELATIVE ERROR 

NUMBER OF ITERATIONS, 
1000 



MEAN 

= 164.285 

= 165.526 

= 138.519 



STAND DEV 
. 1207 
. 1290 
.1075 



BUOY 
. 0080 



TRANS 
. 0080 



TARGET 
. 0080 



TRUE ERROR FOR TARGET 1 
TRUE ERROR FOR TARGET 2 
RELATIVE ERROR 

NUMBER OF ITERATIONS, 
1000 



MEAN 

= 183.286 

= 179.448 

= 154.021 



STAND DEV 
. 1335 
. 1264 
. 1197 



BUOY 

.0090 



TRANS 

.0090 



TARGET 

,0090 



TRUE ERROR FOR TARGET 1 
TRUE ERROR FOR TARGET 2 
RELATIVE ERROR 

NUMBER OF ITERATIONS, 
1000 



MEAN 

= 198.663 

= 199.931 

= 172.593 

BUOY TRANS 

.0100 .0100 



STAND DEV 
.1406 
.1478 
. 1288 

TARGET 

.0100 
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Case number 2 



BUOYLOC . DAT 



0.0, 


0.0, 


5.0, 


0.0, 


0.0, 


1000.0 


3000.0, 


5196.0, 


5.0, 


3000.0, 


5196.0, 


1000.0 


6000.0, 


0.0, 


5.0, 


6000.0, 


0.0, 


1000.0 


9000.0, 


5196.0, 


5.0, 


9000.0, 


5196.0, 


1000.0 



TARGET . DAT 

1000.0, 2000.0, 400.0, 

3000.0, 4236.0, 400.0, 

OUTPUT . DAT 

TRUE ERROR FOR TARGET 1 
TRUE ERROR FOR TARGET 2 
RELATIVE ERROR 

NUMBER OF ITERATIONS, 
1000 



MEAN STAND DEV 

= .001 .0000 

= .002 .0000 

= .001 .0000 

BUOY TRANS TARGET 

.0000 .0000 .0000 



TRUE ERROR FOR TARGET 1 
TRUE ERROR FOR TARGET 2 
RELATIVE ERROR 

NUMBER OF ITERATIONS, 
1000 



MEAN 

23.314 

12.332 

19.779 



STAND DEV 
.0147 
.0058 
.0111 



BUOY TRANS TARGET 

.0010 .0010 .0010 



TRUE ERROR FOR TARGET 1 
TRUE ERROR FOR TARGET 2 
RELATIVE ERROR 

NUMBER OF ITERATIONS, 
1000 



MEAN 

46.881 

24.694 

39.783 



STAND DEV 
.0300 
.0116 
.0226 



BUOY TRANS TARGET 

.0020 .0020 .0020 



TRUE ERROR FOR TARGET 1 
TRUE ERROR FOR TARGET 2 
RELATIVE ERROR 

NUMBER OF ITERATIONS, 
1000 



MEAN 

71.010 

37.106 

60.300 



STAND DEV 
.0467 
.0175 
.0355 



BUOY 

.0030 



TRANS 
. 0030 



TARGET 
. 0030 



TRUE ERROR FOR TARGET 1 
TRUE ERROR FOR TARGET 2 
RELATIVE ERROR 

NUMBER OF ITERATIONS, 
1000 



MEAN 

96.323 

49.592 

81.926 



STAND DEV 
.0673 
.0236 
.0521 



BUOY 

.0040 



TRANS 

.0040 



TARGET 
. 0040 
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TRUE ERROR FOR TARGET 1 
TRUE ERROR FOR TARGET 2 
RELATIVE ERROR 

NUMBER OF ITERATIONS, 
1000 



MEAN 

= 122.966 
= 62.178 
= 104.788 



STAND DEV 
.0913 
.0299 
.0721 



BUOY TRANS TARGET 
.0050 .0050 .0050 



TRUE ERROR FOR TARGET 1 
TRUE ERROR FOR TARGET 2 
RELATIVE ERROR 

NUMBER OF ITERATIONS, 
1000 



MEAN 
148 . 598 
74.547 
127.339 



STAND DEV 
. 1136 
.0364 
.0906 



BUOY 

.0060 



TRANS 

.0060 



TARGET 

.0060 



TRUE ERROR FOR TARGET 1 
TRUE ERROR FOR TARGET 2 
RELATIVE ERROR 

NUMBER OF ITERATIONS, 
1000 



MEAN 

= 174.871 
= 87.047 
= 150.839 



STAND DEV 
.1362 
.0428 
. 1110 



BUOY TRANS TARGET 
.0070 .0070 .0070 



TRUE ERROR FOR TARGET 1 
TRUE ERROR FOR TARGET 2 
RELATIVE ERROR 

NUMBER OF ITERATIONS, 
1000 



MEAN 

= 194.483 
= 99.075 
= 169.387 



STAND DEV 
. 1489 
.0498 
.1218 



BUOY TRANS TARGET- 
.0080 .0080 .0080 



TRUE ERROR FOR TARGET 1 
TRUE ERROR FOR TARGET 2 
RELATIVE ERROR 

NUMBER OF ITERATIONS, 
1000 



MEAN 

= 209.383 

= 111.316 

= 183.449 



STAND DEV 
.1497 
.0572 
.1207 



BUOY TRANS TARGET 
.0090 .0090 .0090 



TRUE ERROR FOR TARGET 1 
TRUE ERROR FOR TARGET 2 
RELATIVE ERROR 

NUMBER OF ITERATIONS, 
1000 



MEAN 

= 232.755 

= 123.187 

= 206.640 



STAND DEV 
.1665 
.0622 
.1412 



BUOY TRANS TARGET 

.0100 .0100 .0100 
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. Case number 3 



BUOYLOC. 


DAT 




0.0, 


0.0, 


5.0 


3000.0, 


5196.0, 


5.0 


6000.0, 


0.0, 


5.0 


9000.0, 


5196.0, 


5.0 



TARGET . DAT 

2036.0, 2036.0, 400.0, 

1965.0, 1965.0, 400.0, 

OUTPUT . DAT 

TRUE ERROR FOR TARGET 1 
TRUE ERROR FOR TARGET 2 
RELATIVE ERROR 

NUMBER OF ITERATIONS, 
1000 



0.0, 


0.0, 


1000 


3000.0, 


5196.0, 


900 


6000.0, 


0.0, 


800 


9000.0, 


5196.0, 


700 



MEAN 


STAND DEV 


. 002 


.0000 


. 001 


.0000 


.001 


.0000 



BUOY TRANS TARGET 
.0000 .0000 .0000 



TRUE ERROR FOR TARGET 1 
TRUE ERROR FOR TARGET 2 
RELATIVE ERROR 

NUMBER OF ITERATIONS, 
1000 



MEAN 

21.890 

21.668 

16.342 



STAND DEV 
.0139 
.0138 
.0095 



BUOY 
. 0010 



TRANS 

.0010 



TARGET 

.0010 



TRUE ERROR FOR TARGET 1 
TRUE ERROR FOR TARGET 2 
RELATIVE ERROR 

NUMBER OF ITERATIONS, 
1000 



MEAN 

44.028 

43.639 

32.941 



STAND DEV 
. 0283 
.0283 
.0193 



BUOY 
. 0020 



TRANS 

.0020 



TARGET 

.0020 



TRUE ERROR FOR TARGET 1 
TRUE ERROR FOR TARGET 2 
RELATIVE ERROR 

NUMBER OF ITERATIONS, 
1000 



MEAN 



STAND DEV 



= 66.771 .0445 

= 66.268 .0446 

= 50.214 .0303 

BUOY TRANS TARGET 
.0030 .0030 .0030 



TRUE ERROR FOR TARGET 1 
TRUE ERROR FOR TARGET 2 
RELATIVE ERROR 

NUMBER OF ITERATIONS, 
1000 



MEAN 

= 90.263 

= 89.303 

= 68.667 

BUOY TRANS 
.0040 .0040 



STAND DEV 
.0629 
.0620 
.0436 

TARGET 
.0040 
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TRUE ERROR FOR TARGET 1 
TRUE ERROR FOR TARGET 2 
RELATIVE ERROR 

NUMBER OF ITERATIONS, 
1000 



MEAN 

= 114.024 
= 114.166 
= 89.502 



STAND DEV 
.0834 
.0862 
. 0640 



BUOY 

.0050 



TRANS 
. 0050 



TARGET 

.0050 



TRUE ERROR FOR TARGET 1 
TRUE ERROR FOR TARGET 2 
RELATIVE ERROR 

NUMBER OF ITERATIONS, 
1000 



MEAN 



STAND DEV 



= 137.825 .1066 

= 136.027 .1048 

= 111.802 .0864 

BUOY TRANS TARGET 
.0060 .0060 .0060 



TRUE ERROR FOR TARGET 1 
TRUE ERROR FOR TARGET 2 
RELATIVE ERROR 

NUMBER OF ITERATIONS, 
1000 



MEAN 

= 164.585 

= 159.147 

= 136.492 

BUOY TRANS 
.0070 .0070 



STAND DEV 
. 1340 
. 1289 
.1275 

TARGET 
.0070 



TRUE ERROR FOR TARGET 1 
TRUE ERROR FOR TARGET 2 
RELATIVE ERROR 

NUMBER OF ITERATIONS, 
1000 



MEAN 



STAND DEV 



= 178.254 .1318 

= 174.554 .1306 

= 146.033 .1085 

BUOY TRANS TARGET 
.0080 .0080 .0080 



TRUE ERROR FOR TARGET 1 
TRUE ERROR FOR TARGET 2 
RELATIVE ERROR 

NUMBER OF ITERATIONS, 
1000 



MEAN 


STAND DEV 


= 201.941 


. 1519 


= 198.101 


. 1477 


= 178.767 


. 1563 


BUOY TRANS TARGET 



.0090 .0090 .0090 



TRUE ERROR FOR TARGET 1 
TRUE ERROR FOR TARGET 2 
RELATIVE ERROR 

NUMBER OF ITERATIONS, 
1000 



MEAN 

= 217.521 

= 212.118 
= 186.845 



STAND DEV 
. 1554 
.1509 
. 1463 



BUOY 
. 0100 



TRANS 

.0100 



TARGET 

.0100 
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Case number 4 



BUOYLOC . DAT 



0.0, 


0.0, 


5.0, 


0.0, 


0.0, 


1000.0 


3000.0, 


5196.0, 


5.0, 


3000.0, 


5196.0, 


900.0 


6000.0, 


0.0, 


5.0, 


6000.0, 


0.0, 


800.0 


9000.0, 


5196.0, 


5.0, 


9000.0, 


5196.0, 


700.0 



TARGET.DAT 

1000.0, 2000.0, 400.0, 

3000.0, 4236.0, 400.0, 

OUTPUT . DAT 

TRUE ERROR FOR TARGET 1 
TRUE ERROR FOR TARGET 2 
RELATIVE ERROR 

NUMBER OF ITERATIONS, 
1000 



MEAN STAND DEV 

= .001 .0000 

= .001 .0000 

= .001 .0000 

BUOY TRANS TARGET 

.0000 .0000 .0000 



TRUE ERROR FOR TARGET 1 
TRUE ERROR FOR TARGET 2 
RELATIVE ERROR 

NUMBER OF ITERATIONS, 
1000 



MEAN 

= 24.415 

= 13.655 

= 20.218 
BUOY TRANS 
.0010 .0010 



STAND DEV 
.0156 
.0068 
.0114 

TARGET 
.0010 



TRUE ERROR FOR TARGET 1 
TRUE ERROR FOR TARGET 2 
RELATIVE ERROR 

NUMBER OF ITERATIONS, 
1000 



MEAN 

49.162 

27.372 

40.721 



STAND DEV 
.0319 
.0136 
. 0233 



BUOY 

.0020 



TRANS 

.0020 



TARGET 

.0020 



TRUE ERROR FOR TARGET 1 
TRUE ERROR FOR TARGET 2 
RELATIVE ERROR 

NUMBER OF ITERATIONS, 
1000 



MEAN 
74 . 694 
41.198 
61.921 



STAND DEV 
. 0504 
.0208 
. 0372 



BUOY 
. 0030 



TRANS 
. 0030 



TARGET 
. 0030 



TRUE ERROR FOR TARGET 1 
TRUE ERROR FOR TARGET 2 
RELATIVE ERROR 

NUMBER OF ITERATIONS, 
1000 



MEAN 

= 102.194 
= 55.192 
= 84.955 



STAND DEV 
.0764 
.0283 
. 0585 



BUOY 
. 0040 



TRANS 
. 0040 



TARGET 
. 0040 
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TRUE ERROR FOR TARGET 1 
TRUE ERROR FOR TARGET 2 
RELATIVE ERROR 

NUMBER OF ITERATIONS, 
1000 



MEAN STAND DEV 
= 130.319 .1024 

= 69.366 .0364 

= 109.312 .0842 

BUOY TRANS TARGET 
.0050 .0050 .0050 



TRUE ERROR FOR TARGET 1 
TRUE ERROR FOR TARGET 2 
RELATIVE ERROR 

NUMBER OF ITERATIONS, 
1000 



MEAN 
154.690 
83 . 145 
130.410 



STAND DEV 
. 1222 
.0447 
.0952 



BUOY 
. 0060 



TRANS 
. 0060 



TARGET 
, 0060 



TRUE ERROR FOR TARGET 1 
TRUE ERROR FOR TARGET 2 
RELATIVE ERROR 

NUMBER OF ITERATIONS, 
1000 



MEAN 
183 . 675 
97.235 
156.063 



STAND DEV 
.1501 
.0547 
.1203 



BUOY 
. 0070 



TRANS 
. 0070 



TARGET 
. 0070 



TRUE ERROR FOR TARGET 1 
TRUE ERROR FOR TARGET 2 
RELATIVE ERROR 

NUMBER OF ITERATIONS, 
1000 



MEAN 

= 202.344 

= 112.773 

= 173.353 

BUOY TRANS 
.0080 .0080 



STAND DEV 
. 1559 
.0695 
. 1238 

TARGET 
.0080 



TRUE ERROR FOR TARGET 1 
TRUE ERROR FOR TARGET 2 
RELATIVE ERROR 

NUMBER OF ITERATIONS, 
1000 



MEAN 

221.332 

126.445 

194.134 



STAND DEV 
. 1593 
. 0794 
. 1413 



BUOY 
. 0090 



TRANS 
. 0090 



TARGET 
, 0090 



TRUE ERROR FOR TARGET 1 
TRUE ERROR FOR TARGET 2 
RELATIVE ERROR 

NUMBER OF ITERATIONS, 
1000 



MEAN 

239.836 

140.192 

210.705 



STAND DEV 
.1740 
.0861 
. 1459 



BUOY 

.0100 



TRANS 

.0100 



TARGET 

,0100 
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